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ANOMALOUS  CONVERGENCE  Or  A CONTINUED  TRACTION 
FOR  RATIOS  OF  KUMMER  FUNCTIONS 
Walter  Gautschi 

1.  Introduction.  From  the  differential  equation  satisfied  by  Kummer's 

2 

function  M(a,b;z)  = 1 + § + ■“TT  + • • ■ > Perron  [4,  p.  278]  de- 

velops the  following  continued  fraction. 


zM'(a,b;z)  _ az  (a-fl)z  (a+2)z 

M(a,b;z)  b-z+  b+l-z+  b+2-z+  * ’ 


b t 0, -1, -2, . 


y 


where  M'(a,b;z)  = (d/dz)M(a , b;z)  = (a/b)M(a+l, b+1  ;z).  While  the  con- 


tinued fraction  converges  for  any  complex  z not  a zero  of  M(a,b;z),  the 
convergence  behavior  can  be  extremely  deceptive,  when  |z|  » max(  | a | , |b|  ), 
particularly  if  Re  z > 0.  The  point  is  illustrated  by  concrete  examples  in- 


volving Bessel  and  incomplete  gamma  functions. 

2.  The  phenomenon  of  apparent  convergence  to  the  wrong  limit. 


We  assume,  for  simplicity,  that  a#0  and  b-z  * 0,  -1, -2, . . . . Equation 


(1.1)  can  then  be  written  in  the  form 

,?  b-z  M'(a,b;z)  J_  _^1  _^2_ 

' ' a M(a,b;z)  1+  1+  1+  * ' * * 

where 


(2.  2) 


(a-t-k)z 

(b-z+k-lXb-z+k) 


k = 1,  2,  3, . . . 


Alternatively  (cf.  , e.g.,  Wall  [6,  p.  17 f f] ) , 


(2.  3) 
where 


b-z  M'(a,b;z)  _ y 
a M(a,b;z)  “ Pk  ’ 


p0  " 1 » Pk  = plp2"*pk’  k = 1, 2,  3, . . . , 
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and 


<2.  4) 


p0  = ° ’ 


~%(l+pk-\} 
pk  ' i + ak(i+pk_1: 


, k = 1,2,3,.  ..  . 


The  infinite  series  in  (2.  3)  represents  the  continued  fraction  in  (2.1)  :n  the 
sense  that  the  n-th  partial  sum  of  the  former  is  equal  to  the  n-th  convergent 
of  the  latter,  n = 1, 2,  3, . . . . 

Evidently,  the  terms  p^  in  the  series  of  (2.  3)  decrease  or  increase 
in  absolute  value  according  to  as  |p^|  < 1 or  |p^|  > 1,  respectively.  It  is 
useful,  therefore,  to  examine  the  behavior  of  | p^  | as  a function  of  k . 

Assuming  I P^  j I < 1 » then  Ip^l  < 1 certainly  if  | a ^ | < i . On  the 

other  hand,  by  (2.2),  if  |z|  > |b|+k,  then  |a  | < ( | a | +k)  | z | ( j z | - | bj  -k)  Z, 

and  an  elementary  calculation  shows  the  upper  bound  for  ! a^  | to  be  <{  if 


(2.  5) 


| z | > 2((b|  + 2 1 a | + 3k)  . 


It  follows  that  (2.  5),  together  with  | p^  <1  , implies  | p^  | <1  . Since, 
initially,  p = 0 , we  obtain  by  induction  that  |p  | < 1 for  all  k satisfying 

U I' 

(2.  5). 

If  j z ] is  large,  we  see  that  the  terms  p^  in  (2.3)  must  decrease 
initially,  the  rate  of  decrease  being  greater  the  larger  |z|  . The  continued 
fraction  in  (2.1)  then  gives  the  appearance  of  converging  rapidly  to  a value 
of  the  order  of  magnitude  1 , yielding  for  M'(a,  b;z)/M(a,  b;z)  a value 
approximately  equal  to  -a/z.  This  is  obviously  the  wrong  answer,  if 
Re  z > 0.  Indeed,  from  known  asymptotic  formulas  [5,  Eq.  13.1.4], 


(2.6) 


M'(a,  b;z) 


M(a,  b;z) 


~ 1 as  I z ! -*  <*>  in  Re  z > 0. 


What  is  likely  to  happen,  then,  is  that  the  terms  p^  , after  the  initial 
descent,  begin  to  increase  again,  and  converge  to  zero  only  after  reaching 
some  peak  values  which  are  sufficiently  large  so  as  to  contribute  to  a limit 
consistent  with  (2.6).  It  is  only  during  the  "final  descent"  of  the  terms  p^ 
that  the  correct  limit  will  be  attained  (assuming  no  rounding  errors). 

The  phenomenon  of  apparent  convergence,  while  prevalent  for  Re  z > 0 
and  | z | large,  need  not  occur  if  Re  z < 0,  since  in  this  case  [5,  Eg.  13.1.6] 


(2.7) 


M'(a , b;z) 
M(a,  b;z) 


- — as  z in  Re  z < 0 . 

z 1 1 


Nevertheless,  we  will  see  in  examples  that  the  phenomenon  persists  if 


- < | arg  z | 


< IT 


albeit  in  a weakened  form. 


3.  The  case  of  real  z . It  is  instructive  to  examine  in  more  detail 
the  case  of  real  arguments  z and  real  parameters  a,b  satisfying  0 < a+1 
if  z > 0 , and  0 < a+1  < b if  z < 0 . 

It  will  be  convenient  to  introduce  the  quantities 


(3.1) 


= 1 + Pb 


k ~ * ' rk  ’ 
for  which  the  recursion  (2.4)  gives 


(3.2) 


<r0  = 1 ’ 


1 


°k  1+a,  <r. 


k = 0, 1, 2, . . . , 


k = 1,2,3, 


k k-1 

3. 1.  The  case  z = x > 0,  a+1  > 0.  We  consider  two  subcases, 


(i)  b-x  > 0,  (ii)  b-x  < 0. 

In  case  (i),  it  follows  from  (2.  2)  that  a^  > 0 for  all  k > 1 , hence 
from  (3.2)  that  0 < cr^  < 1 for  all  k > 1 , and  therefore  from  (3.1)  that 
-1  < < 0.  We  see  that  the  terms  p^  ju  (2.  3)  alternate  in  sign  and  de- 

crease monotonlcally  in  modulus.  In  fact,  since  a^~*  0 , hence  a^<r^  j 0 


-3- 


as  It-**1,  we  have  cr^  — 1,  and  so  p^  — 0,  meaning  that  the  :^ern  i 
converges  faster  than  any  geometric  series.  Indeed, 


(3.  3) 


as  is  easily  verified. 


~ - k as  k~*> 


In  case  (ii),  there  exists  a unique  integer  > 1 such  that 

x-b  c k < x-b+1.  Therefore,  a.  < 0,  while  a.  > 0 for  all  k * k„.  It 
0 kQ  k 0 

follows  as  before  that  -1  < p < 0 for  k < k . If  tr  >0  (even  though 

k0 

a.  < 0 )>  then  -1  < p < 0 also  for  all  k > k , and  we  have  the  same 

0 

alternating  and  supergeometric  convergence  behavior  as  in  case  (i).  If, 

however,  <r  <0  (which  will  be  the  case  if  x is  large),  then  |p  | >1. 

o V 

Since  a^  > 0 for  k > k^,  the  inequality  |p^|  >1  will  persist  as  long  as 

a k rema*ns  negative.  Eventually,  however,  <r^  has  to  turn  positive  (the 

series  in  (2.  3)  being  convergent),  and  from  this  point  on,  all  subsequent  <r's 

remain  positive,  hence  the  corresponding  p's  less  than  one  in  modulus. 

Therefore,  if  x is  large,  the  sequence  { | | } initially  decreases,  then 

increases,  and  finally  decreases  to  zero  at  a supergeometric  rate  given  by 

(3.3).  The  "dip-and-peak"  effect  is  more  pronounced,  the  larger  x,  and  is 

what  gives  rise  to  the  phenomenon  of  apparent  convergence. 

3.  2.  The  case  z = -x  < 0,  0 < a+1  < b.  This  time.  a.  < 0 for  all 
— k 

k > 1.  Noting  that  the  function  x(b+x+k-l) ’^(b+x+k) on  0 < x < oo 

i i 

assumes  a unique  maximum  at  (b+k-1)  2 ( b+k ) 3 , we  find  that 


laJ  5 ~r  , k>i, 

(V  b+k-1  + n I b+k  ) 


and  thus,  in  particular, 


-4- 


a+k 

’^k1  — Zb+Zk-l+Z^b+k-l  \l  b+k 
Since  a+1  < b,  it  follows  that 


a+k 

2b+2k-l+2(b+k-l) 


a+k 
b+k-1  • 


ak<0  - 


| a | < l for  all  k > 1 


From  this,  and  (3.  2),  we  deduce  inductively 

2(  k+1 ) 


1 < cr,  < 
k 


k+2 


for  all  k > 1 , 


hence,  in  particular,  0 < < 1 for  all  k > 1.  The  series  in  (2.  3)  is  now 

a series  of  positive  monotonically  decreasing  terms,  and  convergence 
thus  monotone  and,  as  before,  supergeometric. 

4.  Examples. 

4.1.  Bessel  functions.  We  specialize  (1.1)  to  a = v+|,  b=2v+l, 
where  v > and  use  M( v+j,  2v+l;z)  = r(l+v)exp(|zX4Z)"vMiz),  to- 
gether with  the  differential-difference  relation  I1  (z)  = I ,(z)  - — I (z), 

v v — 1 z v 


to  obtain 


, „-l(*z) 

i (U) 


1 | Vl 


4v 


+ 1 


(4.1) 


(v+-j)z  (v+-|)z  (v+|)z 


2v+l-z+  2v+2-z+  2v+3-z+ 


(v+-)z  » 

2v+l-z  . LJ.  ^k 
k=0 


In  Figure  4.1,  the  moduli  of  the  terms,  | | , are  plotted  in  function  of  k , 

for  v=l,  and  z=  rei<|P,  r = 10,  20, 40,  (p=0,  tt/8,  2^/8, . . . ,ir.  The  behavior 
of  |p  | , when  r is  fixed  and  <p  varies  between  0 and  tt/2,  is  almost 
identical  and  is  represented  by  one  curve  in  Figure  4.1.  The  dependence  on 
9 is  shown  only  in  the  case  r = 40,  but  is  analogous  for  the  other  values  of  r. 
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It  is  seen,  typically,  that  the  terms  [p  | decrease  rapidly  at  the 
beginning,  then  bottom  out  and  rise  to  a sharp  peak,  before  finally  converging 
to  zero.  The  dip  of  the  curve,  and  the  upswing  that  follows,  are  quite  sub- 
stantial if  0<<p  <n/z  and  r is  large,  for  reasons  explained  in  Section  3.1, 
case  (ii).  As  <p  increases  from  it/ 2 to  tt  , the  peaking  of  the  curves  gradu- 
ally weakens  and  finally  disappears  when  <p  = it.  a similar  attenuation  takes 
place  upon  increasing  the  value  of  v , as  is  to  be  expected  from  the  discus- 
sion in  Section  3.1,  case  (i). 

The  seriousness  of  the  convergence  anomaly  can  be  seen,  e.g.  , in 
the  case  v - 1,  r = 40,  0<  cp  < tt/2.  If  we  require  ten  decimal  digit  accuracy, 
we  will  attain  it  at  about  k = 15  and  retain  it  through  about  k = 60,  the  partial 
sums  in  (4.1)  all  having  the  same  value  to  ten  decimal  digits  in  the  range 
15  < k < 60.  This  "apparent  limit",  of  course,  is  totally  incorrect,  as  the 
main  contribution  to  the  series  comes  from  the  few  terms  around  k = 100. 

The  situation  is  aggravated  by  the  fact  that  the  numerical  process  of  gener- 
ating the  terms  p^  is  accompanied  by  a substantial  loss  of  accuracy  during 
the  upswing  of  the  curve,  amounting  to  a loss  of  about  16  digits  when  r = 40. 

A further  complication  is  the  apparent  lack  of  warning  signals:  Known  a- 
posteriori  error  estimates  (see,  e.g.  , [3])  either  do  not  apply,  or  seem  to 
apply  only  in  the  region  of  "final  descent”. 

On  the  other  hand,  the  convergence  behavior  of  the  continued  fraction 
in  (4.1)  is  quite  acceptable  when  <p  - it,  i.e.  , z = -x,  x > 0,  in  which  case 
(4.1)  can  be  given  the  form 
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j 4 v , -<1'+7>X  (^+j)y  C'+jlx 

( 4 ■ Z)  7 X — -1  = - — 

2 j <ix|  x J 2v+l+x-  2 v + 2+x-  2i'+3+x- 

v v 

Convergence  is  more  rapid  the  larger  x and/or  v.  The  use  of  this  con- 
tinued fraction,  in  combination  with  Gauss’  continued  fraction,  is  further 
discussed  in  [2]. 


4.2.  Incomplete  gamma  function.  We  have  M(a,a+l;z)  = 

-a 

a(-z)  \(a,-z),  where  -y(a,-)  denotes  the  incomplete  gamma  function. 

Noting  that  M’(a,a+l;z)  = ~~r M(a+1 , a+2;z),  Equation  (1.1)  now  takes 

a +i 

the  form 


(4.  3) 


\(a+l,  -z)  _ -az  ( a 4-1  )z  (a+2)z 

\(a, -z)  a+l-z+  a+2-z+  a+3-z+  ' 


The  convergence  behavior  of  (4.  3)  for  a > 0 appears  to  be  quite  analogous 
to  that  of  (4.  1),  exhibiting  the  phenomenon  of  apparent  convergence  for  z 
large  in  the  complex  plane  cut  along  the  negative  real  axis.  Along  the 
negative  real  axis,  we  have  monotone  convergence,  if  a > -1,  according 
to  the  results  of  Section  3.  2.  In  this  case, 


(4.  4) 


\(a+l , x) ax  (a+l)x  (a+2)x 

\(a,x) 


a+l+x-  a+2+x-  a + 3+x-  ’ ‘ > 


x > 0. 


We  may  combine  this  with  \(a+i,x)  = a\(a,x)  - x3e'X  to  obtain 


(4.  5) 


-a_x  1 ax  (a+l)x  (at2)x 

x e \(a,x)  = — s— 1 1 — ...  , x > 0 . 

a-  a+l+x-  a+2+x-  a+3+x- 


The  use  of  this  continued  fraction,  in  combination  with  other  methods,  to 
evaluate  incomplete  gamma  functions  is  discussed  in  [1). 
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